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Abstract 

Solution of ground enhancement material application for 
decreasing resistance of grounding system is a current 
problem that draws much attention from scientists 
[(Lightning Protection International PTY), (Van Dinh Anh 
and his Group, 2004)], especially, a grounding system is 
installed in small area or high resistivity of Soil. In order to 
determine the grounding resistances of simple forms 
(horizontal rod or vertical rod. we can use the given 
software (ERICO Company LTD) or the calculating formulas 
(Chuong Ho Van Nhat, 2010). This paper presents some 
other formulas for calculating a grounding resistance of 
simple forms with ground enhancement material. Calculated 
results from the formulas of this paper were verified by the 
results of (Chuong Ho Van Nhat, 2010) and the results from 
the grounding calculating software of ERICO Company in 
Australia. 
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Calculating method 

Grounding Resistance of Single Vertical Rod 

1) GEM Layer of Cylinder Form 

Single vertical rod with diameter d , length / and 
resistivity p3 is buried in the ground with depth t . 
It is environed by the GEM layer with resistivity p2 
and the soil environment with resistivity pi (see 
FIG .1). 

Choosing coordinates system as in FIG.2 and 
considering the potential at point (x,z), according 
to (Nathan Ida, 2000), we have 



<P(z,x) = j" 



Ipdz 



AnUz 1 + x 2 



(1) 



Where / is the current flows through the rod, p is 
the resistivity of environment at point (x, z) and / 
is the length of rod. 




FIG 1 SINGLE VERTICAL ROD IN UNIFORMED ENVIRONMENT 

For simplifying in calculation, the method of 
mirror images is used and the soil is considered as 
a uniformed environment with resistivity pi and 
we do not mention influence of temperature upon 
the environment containing grounding system. So, 
the potential along the rod is calculated as follows: 



<p(x)=\ d<p(z,x)= d<p(z, x)+\ d<p(z,x) 

J—l/2, J J2.t— 1/2 



4x1 
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-Arsh 



2t-U2 



Replacing Arsh (x)=ln [(x)h- ^/jc^+lj and similarly 
for other variables of (2), we receive: 



q(x)-- 



4x1 



21n 



' ^f+Ax 1 +l\ A ^{4t+lf+4x 2 +4t+l 



2x 



^(4t-lf+4x 2 +4t-l 



(3) 



Grounding Resistance of Vertical Rod with the GEM 

This resistance consists of three components (see FIG. 
2): 
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FIG.2 GROUNDING ROD WITH GEM 
1) The Component in the Metal Environment 
The potential function 

Acp l = (p{Q) - (p{d 1 2) 
with p = p3. From (3), we have: 
Ip 3 



and 



<p(d/2) = 



Ad 



IPs 
2M 



In 



r At+0 
At-l 



In 
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+ ii„ 
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From (4) and (5), we receive 



A(p, = 



Ip 3 



2jd 



In 



21 



\U J 



So, the resistance value of this component 

'2/V 



R, 



= A 
/ 2^ 



In 



2) The Component in the GEM Environment 
The potential function 

A<p 2 = <p{d 1 2) - <p{C) 

with p = p2. 

when x= d 12, x=C and note that 



« (At + l) 2 , we have: 



V ■ 



a 7 A 



In 



d 2 Uf-/ 



V(/ 2 +4C 2 )+/^|_ 1 J ^{At+lf +AC 2 +4t+l 



2C 
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and 



A<p. 



(4) 



(5) 



(6) 



(7) 
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3) The Component in the Soil Environment 

The potential function A<p 3 = <p(C) - <p(x>) with p = 
pi. Because of 9(00) = 0, A<p 3 = (p(C) ■ So: 



21r 



2C 



+lr 



>(At+lf+AC ! +At+l 



A^ 3 



2tf 
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2 
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(10) 



(11) 



At last, the total resistance value: 
R C =R I + R 2 + R 3 

In reality, a thickness of the chemical layer C is 
very smaller than the rod length / (C« / ). 
Therefore, from (7), (9) and (11), we have: 

f IsjAt + I A 



R, = 



2 nl 



Mi* 



C-jAt-l 



(12) 



Grounding Resistance of Equivalent Vertical Rod 

For calculating grounding resistance of simple forms 
with or without the GEM, we wish to use only a form 
of a common formula. So, the formula for calculating 
such equivalent grounding resistance value will be 
found in this paper. 

We can consider that the single vertical rod is only 
buried in one single soil layer consisting of GEM and 
natural soil with their equivalent resistivity Q eq . This 
resistance value consists of only two components (see 
FIG.3). 

Calculating similarly as the above points, we have: 
1) The Component in the Metal Environment 



*i'= 



2nl 



In 



(13) 



2) The Component in the Equivalent Soil 
Environment 



Peg 

Anl 



21n| ^ |h In 



r (At + Z)^ 
(4?-/) 



(14) 
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FIG .3 VERTICAL GROUNDING ROD AND ITS IMAGE 
Total resistance: Rc = R'i + R'3 
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Hence, In order to determine the equivalent 
resistivity we will equalize (12) and (15), from this 
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The equivalent soil resistivity will be 
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(16) 



Hence, we can count as single grounding vertical 
rod (with parameters d ,l,p3) that is only buried in 
one single soil layer with equivalent resistivity Q eq . 
This equivalent resistivity is calculated by equation 
(16) (see FIG.4). 



fmrnmnmwm 




i 

FIG .4 VERTICAL ROD 
IN EQUIVALENT SOIL ENVIRONMENT 

In reality, p 3 « 0, we have from (15): 

P eq 



R, 



2nl 



r 21 1 , 

In — + —In 

d 2 



f (4f + I) A 



(4r - Z) 



(17) 



with p g? is calculated by (16) . If there is not the 

GEM layer, then p = p x and (17) will become: 

(see (Alya B. Joffe, 2010), (ANSI/IEEE Std 80-2000), 
(Chuong Ho Van Nhat, 2009), (Viet Hoang, 2010)) 



R, 



2 Til 



21 1, 
In — + — In 

d 2 



f {At + I) ^ 
(At - I) 



(18) 



3) GEM Layer of Rectangular Prism Form 

We consider the single vertical rod with diameter 

^ , length ^ and resistivity p3. It is surrounded by 
the GEM layer with resistivity p2 and the soil 
environment with resistivity pi. It is fixed in the 
ground with the depth t and the GEM layer has the 

dimensions ax bxl ( see FIG5) 




FIG.5 SINGLE GROUNDING VERTICAL ROD WITH THE GEM 
LAYER 

In this case, we can convert the rectangular prism 
form into a cylinder form with its diameter C [7] 
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(19) 



and we will receive the formulas as in the case of 
the GEM layer of cylinder form. 

Grounding Resistance of Single Horizontal Rod 

1) GEM Layer of Cylinder Form 

Single horizontal rod with diameter d , length / 
and resistivity p3 is fixed in the ground with the 
depth t . It is environed by the GEM environment 
with resistivity pi and the soil environment with 
resistivity pi (see FIG. 6) 




FIG.6 SINGLE GROUNDING HORIZONTAL ROD WITH THE 
GEM LAYER 

Applying (1) and after some transformations, we 
found the magnetic field that is created by the 
horizontal rod at the point (x, z) 



<Pi(x) = T^rln 
2nl 



2x 



2x 



+ 1 



(20) 



and the magnetic field is created by the image of 
the horizontal rod at the point (x, z): 



cp„ (x) = — — In 
2nl 



I 



2(2t-x)x 



I 



2(2t-x) 



(21) 



Grounding Resistance of Horizontal Rod with the 
GEM 

We use the calculating method as at the above 
point and receive the results as follows : 

1) The Component in the Metal Environment 

A <P> _ Ps J dT\ (22) 



(23) 



Ri = 

2) The Component in the GEM Environment 



A< Pi _ Pi ln f C(2t - C) 



2nl 



dt 



3) The Component in the Soil Environment 



_ IPx 

I 2nl 



In 



C(2t-C) 



(24) 



So, total resistance: Rt = Ri + R2 +R3 



Rj — 



2nl 



(A)lnf— — 
V " \C(2t-C) 



(25) 



f C(2t-Q ] , , d 
+ p 2 H— — +A ln- 



Grounding Resistance of Equivalent Horizontal Rod 

We use the calculating method as at the above 
point and receive the results as follows: 

1) The Component in the Metal Environment 



._ Hi 



2nl 



In 



21 



(26) 



2) The Component in the Equivalent Soil 
Environment 



In 



2nl 

Total resistance: Rt = R'i + R'3 
1 



y dt j 



Rj — 



2nl 



, I d 

, + p 3 In — 

dt I 3 {21 



(27) 



(28) 



Hence, In order to determine the equivalent 
resistivity we will equalize (25) and (28), from this 



A In 



1- 



C(2t-C) 



+ p 2 ln 



C(2t-C 
dt 

~ P eq In 
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The equivalent soil resistivity will be 



ln 



Pea =" 



I 2 



C(2t-C) 



C(2t-C) 
dt 



(29) 



dt 



So, we can count as single grounding horizontal 
rod (with parameters d , / ,p3) that is only buried in 
one single soil layer with equivalent resistivity Q eq . 
This equivalent resistivity is calculated by equation 
(29) (see FIG.7). 

In reality, p 3 ~ , from (28) we have: 



P eq 



Inl 



ln 



K dt j 



(30) 
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FIG.7 GROUNDING HORIZONTAL ROD IN EQUIVALENT 
SOIL ENVIRONMENT 

When there is not GEM, then p = p l and the 
formula (30) will become: (see [1], [2], [3] and [10]) 



R 



Pi 



2kI 



In 



v dt j 



(31) 



3) GEM Layer of Rectangular Prism Form 




FIG.8 SINGLE HORIZONTAL ROD WITH THE GEM 

Similar as above, we can convert a rod of the 
rectangular prism form (see FIG.8) into a cylinder 
with diameter C is calculated by formula (19) and 
receive the formulas as in the case of the cylinder 
form GEM layer. 

In summary, in order to calculate a grounding 
resistance of simple vertical rod, we use the 
formula 



R, 



2 Til 

without GEM then 
with GEM then 



.21 1 , 
In — + — In 

d 2 



r (At + 0^ 



(At -I) 



Pea = Pl 



In 



Pea = 



I si At + I 
CsjAt-l 



2C 




and in order to calculate a grounding resistance of 
simple horizontal rod, we use the formula 



Rj — 

2 Til 

without GEM then p e = p 



In 



ydt j 



with GEM then 



In 



Pea =- 



I 2 



C(2t-C) 



( C(2t-C) 
\ dt 



In 



ydtj 



Calculating application 

Calculating Data 

We use the following data: resistivity of soil 200 Qm ; 
resistivity of GEM 20 QCm; resistivity of metal rod 
0.000172 QCm; diameter of GEM layer of cylinder 
form 20Cm ; side of GEM layer of square form 20Cm ; 
the depth of horizontal or vertical rod 50Cm. 

Comparing Results Between This Paper and Formulas 
of (Chuong Ho Van Nhat, 2010) 

1) Surveying Vertical Rod 

TABLE 1 COMPARING RESULTS BETWEEN THIS PAPER AND 
FORMULAS OF (Chuong Ho Van Nhat, 2010) 



Length (m) 


3 


4 


5 


6 


7 


8 


9 


10 


Paper 


47.6 


35.9 


29.3 


25.0 


21.8 


19.4 


17.6 


16.1 


(Chuong Ho 
Van Nhat, 2010) 


47.6 


35.9 


29.3 


25.0 


21.8 


19.4 


17.6 


16.1 



TABLE 2 COMPARING RESULTS BETWEEN THIS PAPER AND 
FORMULAS OF (Chuong Ho Van Nhat, 2010) 



Length(m) 


3 


4 


5 


6 


7 


8 


9 


10 


Paper 


46.4 


35.8 


28.6 


24.3 


21.3 


19.0 


17.2 


15.7 


(Chuong 
Ho Van 
Nhat, 2010) 


46.4 


35.8 


28.6 


24.3 


21.3 


19.0 


17.2 


15.7 



2) Surveying Horizontal Rod 

TABLE 3 COMPARING RESULTS BETWEEN THIS PAPER AND 
FORMULAS OF (Chuong Ho Van Nhat, 2010) 



Length(m) 


5 


10 


15 


20 


25 


30 


35 


40 


Paper 


35.8 


22.3 


16.6 


13.3 


11.2 


9.7 


8.6 


7.7 


(Chuong 
Ho Van 
Nhat, 2010) 


35.8 


22.3 


16.6 


13.3 


11.2 


9.7 


8.6 


7.7 
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TABLE 4. COMPARING RESULTS BETWEEN THIS PAPER AND 
FORMULAS OF (Chuong Ho Van Nhat, 2010) 



Length(m) 


5 


10 


15 


20 


25 


30 


35 


40 


Paper 


35.1 


21.9 


16.3 


13.2 


11.1 


9.6 


8.5 


7.7 


(Chuong Ho 
Van Nhat, 
2010) 


35.1 


21.9 


16.3 


13.2 


11.1 


9.6 


8.5 


7.7 



3) Remark 

From the data in TABLES 1, 2, 3 and 4 we realize 
that the formulas of this paper and (Chuong Ho Van 
Nhat, 2010) give the same results. 

Comparing Results between This Paper and GEM 
Software ofERCO Company 

1) Surveying Vertical Rod 

TABLE 5 COMPARING RESULTS BETWEEN THIS PAPER AND 
SOFTWARE OF ERICO COMPANY 



Length(m) 


3 


4 


5 


6 


7 


8 


9 


10 


Paper 


47.6 


35.9 


29.3 


25.0 


21.8 


19.4 


17.6 


16.1 


Software 


40.2 


32.5 


27.4 


23.8 


21.1 


19.0 


17.3 


15.9 


TABLE 6 COMPARING RESULTS BETWEEN THIS PAPER AND 
SOFTWARE OF ERICO COMPANY 


Length(m) 


3 


4 


5 


6 


7 


8 


9 


10 


Paper 


46.4 


35.8 


28.6 


23.1 


21.3 


19.0 


17.2 


15.7 


Software 


39.0 


31.5 


26.6 


23.1 


20.6 


18.5 


16.9 


15.5 



2) Surveying Horizontal Rod 



TABLE 7 COMPARING RESULTS BETWEEN THIS PAPER AND 
SOFTWARE OF ERICO COMPANY 



Length(m) 


5 


10 


15 


20 


25 


30 


35 


40 


Paper 


35.8 


22.3 


16.6 


13.3 


11.2 


9.7 


8.6 


7.7 


Software 


32.1 


20.4 


15.4 


12.4 


10.5 


9.1 


8.1 


7.3 


TABLE 8. COMPARING RESULTS BETWEEN THIS PAPER AND 
SOFTWARE OF ERICO COMPANY 


Length(m) 


5 


10 


15 


20 


25 


30 


35 


40 


Paper 


35.1 


21.9 


16.3 


13.2 


11.1 


9.6 


8.5 


7.7 


Software 


31.3 


20.1 


15.1 


12.2 


10.4 


9.0 


8.0 


7.2 



3) Remark 

From the data in TABLES 5, 6, 7 and 8 we realize 
that the formulas of this paper and the software of 
ERICO Company give the calculating results with 
the errors that are not great (about 15%) with short 
length of rod and are approximate with long. 

Conclusions 

1) . We received other formulas for calculating the 
resistance of simple grounding forms with influence of 
the GEM in comparison with (Chuong Ho Van Nhat, 
2010). 

2) . We received formulas for calculating the equivalent 
resistivity of the GEM and natural soil layer. 

3) . We can use a common formula for calculating the 
resistance of simple grounding forms with or without 
the GEM. 

REFERENCES 

Alya B. Joffe, Kai-sang Lock. "Grounds for Grounding: A 

Circuit to System Handbook". Wiley. -2010. 
ANSI/IEEE Std 80-2000. "IEEE guide for safety in AC 

substation grounding". New York. 
Chuong Ho Van Nhat. "Problems in High Voltage 

Engineering". Publisher HoChiMinh City National 

University, VietNam, 2009. 
Chuong Ho Van Nhat. "Calculating resistance of simple 

grounding forms with or without the soil improved 

chemical substance". 2010 Asia Pacific Power and Energy 

Engineering Conference (APPEEC), March 28-31, 2010. 

Chengdu, China. 
ERICO Company LTD. "The GEM software for calculating 

grounding resistance". 
Lightning Protection International PTY.- LTD. Global 

Lightning Technologies.- Hobrut, Tasmania, Australia 

7000. 

M. E. Ilrusalimov; N. N. Orlov. "High Voltage Engineering". 

Publisher Ralianskaia. Kiev, Ucraine. 
Nathan Ida. "Engineering Electromagnetics". Hamilton 

Printing Co., Rensselaser, NY., 2000. 
Van Dinh Anh and his Group. "Researching and suggesting 

some solutions for improving grounding system of the 



30 



International Journal of Automation and Power Engineering (IJAPE) Volume 2 Issue 2, February 2013 



www.ijape.org 



transmission line and high voltage transformer 
substation of HoChiMinh City Power Company". 2004. 
Viet Hoang. "High Voltage Engineering: Overvoltages in 
Power System". Publisher HoChiMinh City National 
University. VietNam, 2010. 

Authors Introduction 

Chuong Ho Van Nhat was born in 
1954, Hue province, Vietnam. He 
received bachelor of science degree in 
electrical engineering from Ho Chi 
Minh City University of Technology 
in 1977 and Philosophy Doctor degree 
in High Voltage Engineering from 
Kiev Polytechnique University, 
Ukraine in 1997. Now, he has been as 
senior lecturer of Power Systems Department of HCMCity 
University of Technology. His main research interests 
concentrate on the fields of Grounding System, Overvoltages 
in Power Systems and Measurements on High Voltage. 




Tung Nguyen Thanh was born in 1980, 
j^B^k TraVinh province, Vietmam. He received 
bachelor of science degree in electrical 
and electronics engineering in 2007 and 
'^^ta master °f science degree in 2012 from 
H fl University of Technical Education 
ft j M HCMC, HCM City, Vietnam in 2007. He 
H^rfl I has been as lecturer of Dong Thap 
Vocational College. His main research interests concentrate 
on the fields of Grounding System and Renewal Energy. 

Lan Ngo Kim was born in 1976, Vinh 
Phuc province, Vietnam. He received 
bachelor of science degree in electrical and 
electronics engineering from Ho Chi Minh 
City University of Technology - HUTECH 
in 2001 and and master of science degree 
in electrical engineering from University 
of Technical Education HCMC, HCM 
City, Vietnam in 2009. He has been as lecturer of Dong Nai 
Vocational College. His main research interests concentrate 
on the fields of Grounding and Automation in Power 
Systems. 




31 



